Static and dynamic properties of a reversible gel 



P.I. Hurtado, 1 P. Chaudhuri, 2 L. Berthier, 2 and W. Kob 2 

1 Institute Carlos I for Theoretical and Computational Physics, and Departamento de Electromagnetismo y Fi'sica de la 
Materia, Universidad de Granada, Granada 18071, Spain 

2 Laboratoire des Coloides, Verres et Nanomateriaux, Universite Montpellier 2 and CNRS, Montpellier 34095, 
France 

We study a microscopically realistic model of a physical gel and use computer simulations to 
investigate its static and dynamic properties at thermal equilibrium. The phase diagram comprises 
a sol phase, a coexistence region ending at a critical point, a gelation line, and an equilibrium 
00 ' gel phase unrelated to phase separation. The global structure of the gel is homogeneous, but 

C _ 1 the stress is supported by a fractal network. Gelation results in a dramatic slowing down of 

' the dynamics, which can be used to locate the transition, which otherwise shows no structural 

£\J i signatures. Moreover, the equilibrium gel dynamics is highly heterogeneous as a result of the 

presence of particle families with different mobilities. An analysis of gel dynamics in terms of 
^ ■ mobile and arrested particles allows us to elucidate several differences between the dynamics of 

O ' equilibrium gels and that of glass-formers. 
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I. INTRODUCTION 



A gel can be roughly defined as a low-density, disordered, solid material composed by a liquid matrix in which 
dispersed particles form a very open network. In this way, a gel can be notably elastic and jellylike, as for instance 
gelatin, or rather solid and rigid, as silica gel (1). Though gels are materials common to everyday experience, their 
structural and dynamical properties remain puzzling in many respects (2). This is mostly due to the wide window of 
timescales and lengthscales which determine their physical behavior, e.g. from the molecular size of particles in the 
solvent to supramolecular structures (1). An important parameter for a gel is the typical lifetime of the interparticle 
I ■ bonds which define the underlying stress-sustaining network. In the limit of very strong, permanent bonds (with 
typical energies much larger than fc^T), diffusion- limited cluster aggregation (DLCA) may lead to the formation of a 
fractal, system-spanning network usually termed as chemical gel and whose properties follow directly from geometry. 
In particular, gelation in chemical gels can be unambiguously identified with percolation. On the other hand, weaker 
bonds (with energies competing directly with fcgT) result in the formation of physical (or reversible) gels, where the 
links have a finite lifetime and the transient character of the network results in a complex interplay between structure 
and dynamics, leading to non-trivial flow properties. 

Despite its ubiquity, the nature of physical gelation is still under intense debate and several mechanisms have been 
proposed to account for this transition, among them geometric percolation (3-5), the glass transition (6-8), and 
■^j- . arrested phase separation (9; 10). For instance, detailed experiments performed with colloidal particles with tunable 
interactions (11) revealed that a non-trivial interplay between phase separation and kinetic arrest may produce gel-like 
structures. Associating polymers constitute another well-studied example of reversible gels (1). In that case, gels can 
be obtained far from phase separation, producing viscoelastic materials with highly non-linear rheological properties 
that are not well understood (12; 13). In fact, one may speak about different routes to gelation (2): a nonequilibrium 
■ route based on kinetic arrest during spinodal decomposition (an irreversible process), and an equilibrium route in 
which the gel state is reached from an ergodic phase always allowing for the equilibration of the underlying network 
structure (as in the case of polymer association). In all cases, the emergence of the gel phase ensues some degree 
of kinetic arrest (at least for low wave-vector modes) which reflects the appearance of a percolating macroscopic 
structure capable of slowing down particle motions over long timescales. This kinetic arrest implies in many cases a 
close similarity between gelation and glass formation (1). We explain below this similarity for the case of equilibrium 
gels, but discuss also important differences. 

In this paper we report results for a recently introduced model of a reversible physical gel which is microscopically 
realistic (we are in fact inspired by one particular material), complementing the results already published in Ref. 
(3). Moreover, we describe in detail a hybrid Monte Carlo / molecular dynamics numerical approach designed to 
successfully bridge the gap between microscopic details and macroscopic observations, while offering deep insight on 
the nature of physical gels. 



II. A MODEL OF REVERSIBLE GELATION 



Our model is inspired by a specific material (13): A microemulsion of stable and monodisperse oil droplets in water 
mixed with telechelic polymers. These polymers are long hydrophilic chains with hydrophobic end caps. Hydropho- 
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FIG. 1 Left: Sketch of the oil-in-water micro-emulsion and our modelling strategy. Both bridging polymers and the coarse- 
grained effective links used to model them are depicted. The thicker and dashed links represent a Monte Carlo move for the 
connectivity matrix. Right: Potential energy for two droplets connected by a polymer, for different amplitudes 62. In this 
paper we choose £2 = 50 (in blue). Notice the short-range repulsion at r = a — 1, and the divergence at the maximum polymer 
extension, I = 3.5a. 



bicity guarantees that a polymer can form a loop around a single droplet or, more interestingly, a bridge between two 
droplets, see the sketch in Fig. 1 (left). Polymer bridges induce an effective entropic attraction between the droplets 
they connect, due to the confrgurational constraint resulting from the two polymer end caps lying at a given distance. 
For sufficiently high polymer and droplet concentrations, a network of connected droplets spanning the entire system 
can be formed and the system becomes a soft solid. However, thermally activated extraction of the hydrophobic 
heads leads to a slow reorganization of the network structure, and the material eventually flows at long times: it is a 
transient network fluid (14). 

This material is interesting because functionality, lifetime of the bonds, volume fraction, strength of the links can 
all be adjusted independently, which is not always possible in attractive colloids (11), or in previous model systems, 
unless specific adhoc assumptions are made (6-9). For instance, attractive colloidal gels can exhibit long- living bonds 
whenever a large attraction strength between the particles is present, but such a strong attraction simultaneously 
leads to liquid-gas phase separation preventing the formation of a gel at thermal equilibrium (11). On the other hand, 
in our model system both the bond lifetime and the attraction strength can be independently adjusted, allowing a 
different route to colloidal gelation, as discussed below. 

Modelling such a complicated self-assembly of solvent, droplets and complex polymers is a challenge because of 
the wide range of scales involved. In our simplified model we neglect the solvent (which might only affect short time 
dynamics) and include polymers and droplets as the elementary objects. Moreover, since the internal dynamics of 
the polymers (~ 10 -8 s) is much faster than the gel dynamics (~ 1 s), we coarse-grain the polymer description and 
only retain their effect as links inducing an effective interaction between the two droplets they connect, see left panel 
in Fig. 1. Coarse-graining is a crucial step for efficient large scale simulations, not used in previous models (15). 

We consider an assembly of N droplets of diameter a and mass m, interacting, in the absence of polymers, via a 
pair potential typical of soft spheres 

where is the distance between droplets i and j, and t\ is an energy scale. The potential is cut off and regularized 
at a finite distance 2.5 a. In addition, N p polymers of maximal extension I can form bridges between droplets, or 
loops. We assume that polymer loops have an energy cost eo, but no effect on droplet dynamics. On the other hand, 
bridging polymers induce an entropic attraction between connected droplets, which we model using the classic FENE 
form, which is obtained assuming a simple random walk description for polymers but retains the idea that polymers 
have a maximal possible extension (1): 

ra<t, (1) 

so that polymers act as linear springs at small elongation, but cannot become longer than £. The right panel in Fig. 1 
shows the total interaction potential for two droplets of diameter a = 1, linked by one telechelic polymer at different 
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FIG. 2 Phase diagram for a wide range of volume fractions, <j>, and number of polymer heads per droplet, R. Symbols are 
the investigated state points in the sol (O)j S e l (P), and phase separated (o) regions. A yellow point shows the approximate 
location of the critical point. Transition lines are sketched. 

values of €2- A system configuration is thus specified by the droplets positions and velocities, {rj(i), v.;(i)}, and by 
the polymer N x N connectivity matrix, {Cy}, where Cy is the number of polymers connecting droplets i and j. 
The Hamiltonian is thus 

N 

n = E (? v " + + E ^ ) + c ^ ^ ^ )] ) ■ ( 2 ) 

i=l ~ j>i 

Simulation proceeds by solving Newton's equations for the droplets. Lengthscalcs arc given in units of a, energy 
in units of e\, and times in units of \J ma 2 /e\. We use the velocity Verlet algorithm with discretization h — 0.005. 
Simultaneously, we use Monte Carlo dynamics to evolve the polymers. In one elementary move, a polymer is chosen at 
random, and one of its end caps is moved to a randomly chosen neighboring droplet. This proposed move is accepted 
with a Metropolis rate 

r i7nk min 

where AV2 is the potential energy change during the move, T is the temperature, fcs the Boltzmann constant, and 
Tiink controls the timescale for polymer rearrangements. This process is sketched in Fig. 1 (left). Notice that the 
droplets configuration remains unchanged during a polymer move, an approximation justified by the broad separation 
of time scales between the fast polymer dynamics and the much slower droplet motion. In experiments, Tii n k has 
an Arrhenius behavior associated to the excitation cost for polymer extraction. We set I — 3.5 a as measured in 
experiments (13), T = 1, and eo = 1. We found little influence of £2 on the phase diagram, though the choice £2 = 50 
seems most appropriate for a proper numerical comparison with experiments (see Fig. 1, right). The relevant control 
parameters for thermodynamics are the droplet volume fraction, <fi = ira 3 N/(6V), where V is the volume, and the 
number of polymer heads per droplet, R = 2N p /N. Additionally, the dynamics is crucially affected by rii n k, which 
controls the rate for polymer extraction but has no influence on static properties. We performed simulations for a 
wide range of parameters, R £ [0, 18], <j> £ [0.01,0.3], r link E [1, 10 4 ], N E [10 3 , 10 4 ]. 

III. PHASE DIAGRAM AND STRUCTURE 

The phase diagram, as obtained after a systematic exploration of the control parameter space, is shown in Fig. 2. 
Its topology is in good agreement with experiments (13). It comprises a sol phase, a coexistence region ending at a 
critical point, a gelation line determined by geometric percolation, and an equilibrium gel phase. 

In the \ow-cj), low-i? region, the system resembles a dilute assembly of soft spheres, and has the properties of a 
simple liquid. This is the sol phase, see Fig. 3-a. There the largest cluster of droplets connected by polymers (shown 
in light gray in Fig. 3) does not percolate across the system. Increasing R increases the effective attraction between 
droplets, so that phase separation occurs at large R between a \ow-<j>, low-i? phase and a large-0, large- R phase (16), 
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FIG. 3 Simulation snapshots for N = 10 4 . Particles colored in light gray belong to the largest cluster, (a) Sol phase, <j> = 0.15, 
R — 2. (b) Nucleation regime, <f) = 0.01, R = 8. (c) Spinodal decomposition, (j> — 0.05, R = 8. (d) Structure of the physical 
gel for cj> = 0.2, R — 4. (e) Gel right at percolation, cf> = 0.2, R — 2. A percolating (light gray) cluster of droplets connected by 
(red) polymers, through which the remaining (dark blue) droplets can diffuse, see zoom (f). 

see Figs. 3-b and 3-c. We detect phase coexistence from direct inspection of typical configurations in obvious cases. 
We also measured the static structure factor 

/ 1 N N \ 

\ fe=l j=l I 

which exhibits a q~ A scaling behavior at small q usual in coarsening systems, and we also use the so-called demixing 
or inhomogeneous parameter tp n . This is defined by dividing the system in n 3 boxes and measuring the average local 
excess density 

1 " 3 

i>n = — Y] [pi - (p)] , 

1=1 

where pi is the local density in the i-th box, and (p) — N/V = Q<j)/(iicr z ) is the average system density. For a 
homogeneous phase the demixing parameter is zero, but it grows when density inhomogeneities appear, as those char- 
acterizing liquid-gas coexistence. Using these observables we detect both nucleation and spinodal regimes depending 
on the quench depth, see Figs. 3-b and 3-c, respectively. 

For 4> > 0.05 and moderate values of R a broad equilibrium gel region exists between the sol phase at low-i? and 
phase separation at large- R, see right panel in Fig. 2. In the gel phase, a system-spanning homogeneous cluster of 
polymer-connected droplets emerges, Fig. 3-d, which endows the fluid with viscoelastic properties. The observation 
of an equilibrium gel regime in the phase diagram results from the independence of the lifetime and the strength of 
the polymer bonds linking the droplets, and makes this model very interesting to test at once different and usually 
contradictory theoretical frameworks for the gelation transition (2). 

The low-i? sol-gel transition coincides in our model with geometric percolation of polymer-connected clusters. The 
nature of this gelation transition is geometric, rather than thermodynamic, but it cannot be directly detected in 
experiments with structural probes. In fact, the structure of the system does not change across the transition, as 
detected by the structure factor S(q), which remains typical of a simple homogeneous fluid, see Fig. 3-d. This 
is also observed in experiments (13). However, despite this apparent homogeneity, the spanning cluster is highly 
fractal near percolation, and becomes thicker deeper in the gel phase. This can be detected in the low wavevector 
behavior of Si c (q), the structure factor restricted to droplets belonging to the largest cluster, which exhibit power-law 
scaling at low-g with an exponent related to the cluster fractal dimension. On the other hand, at percolation the 
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FIG. 4 Distribution of cluster sizes for N = 10 4 , = 0.2 and increasing values of R across percolation. Notice the power-law 
behavior P(n) ~ n -2 ' 18 near percolation (shown as a dashed line), in agreement with random bond percolation, and the peak 
for large cluster sizes above percolation. 
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FIG. 5 Left: Connectivity distribution Il(fc) for N = 10 4 , (j> = 0.1 and different values of R. Notice that for large R the 
distribution becomes bimodal, with a peak at k = 0, reflecting phase coexistence. The inset shows IT(fc) for R = 4 and 
increasing values of <j>. Right: I4(fc) and its decomposition for particles which belong to the largest cluster or not, for <j> = 0.1 
and different values of R representative of the three phases in Fig. 2. 



distribution of cluster sizes, P(n), shows scaling behavior typical of the random bond percolation universality class. 
In Fig. 4, we show P(n) for = 0.2 and different values of R across percolation. As we approach percolation, the 
cluster distribution develops an algebraic tail with an exponential cutoff which shifts to larger and larger values as the 
transition is approached. Above the transition, a peak for large cluster sizes appears due to the presence of a system- 
spanning macroscopic cluster. Near percolation, we find P(n) ~ n~ 7 , with 7 w 2.18 compatible with random bond 
percolation, as seen in other systems (4; 6; 8; 9). Moreover, the crossover in the scaling of P(n) allows us to determine 
the critical polymer concentration quite accurately. For (j> = 0.2, we find this critical value to be R p (<fr) = 1.85(5). 
Unfortunately, typical experiments have no access to both Si c (q) and P(n) because they cannot discriminate between 
different clusters, so no structural signature of the equilibrium gelation transition is actually detected in experiments. 
However, we show below that the dynamics of the system is strongly affected by percolation, which can then be used 
to locate the transition. 

The coexistence line ends at a second-order critical point, see Fig. 2. The structure of the system becomes highly 
correlated when approaching this liquid-gas phase transition. At this point, we find that S(q) develops an algebraic 
tail with an exponent close to —1.5 at small q. We made no effort in determining accurately the critical point, but 
broadly estimated its location near <f> c w 0.05 and R c « 3.5. 

We have also studied the statistics of droplet connectivity due to polymer bridges. We define II(fc) as the connectivity 
distribution, i.e. the probability density for finding a droplet with k linked neighbors. Figure 5 (left) shows II(fc) 
measured for (j> = 0.1 and different values of R representative of the three phases in Fig. 2. In all cases the connectivity 
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distribution decays exponentially fast for large values of k, with a typical connectivity scale which increases with R. 
One can also detect a crossover from unimodal to bimodal behavior in II(fc) as the coexistence line is crossed. In 
particular, for large R the distribution II(fc) shows a weak peak at k = associated to completely disconnected 
droplets freely diffusing in the low-density phase. Such unimodal-bimodal crossover in II(fc) is a further evidence for 
phase coexistence. The inset to the left panel in Fig. 5 shows the (^-dependence of II(fc) for a fixed value of R = 4. 
Interestingly in the low-<p regime studied the connectivity distribution becomes more homogeneous (i.e. peaked) as 
(f> increases. This might be due to the influence of the liquid-gas critical point. For cf) = 0.1 and R = 4, which is 
relatively close to the critical region, large density fluctuations associated to the critical region are still apparent, and 
this inhomogeneity in the droplet structure is reflected in a more heterogeneous degree distribution (see also below). 
However, increasing <f> moves the system away from the critical region, the system becomes more homogeneous and 
so does the degree distribution. 

We can also decompose the connectivity distribution between a contribution from the largest cluster present in the 
system and its complement. This is done in the right panel in Fig. 5 for <j> = 0.1 and different values of R representative 
of the three phases (sol, gel, phase coexistence). In all cases we detect a clear asymmetry between the connectivity 
distribution of the largest cluster and the rest of the system: most of the highly connected droplets belong to the 
largest cluster. This also means that the largest cluster concentrates more polymer bridges than the average. Such 
asymmetry is most important in the coexistence regime, see Figs. 5 (right) and 3.b-c, where most bridging polymers 
are concentrated in the high-density droplet phase. 

The above discussion concerns relatively low volume fractions, <j> < 0.30. Interesting questions arise when pondering 
over the behavior of our model system at higher values of <fr. At higher volume fraction, we anticipate that the gel 
dynamics described in this work will compete with the kinetic arrest due to the approach of the glass transition 
usually observed in colloidal systems at high density. 

IV. SLOW AND HETEROGENEOUS DYNAMICS 

In this section we study the dynamics of the equilibrium gel, showing that the gel phase indeed behaves dynamically 
as a soft viscoelastic fluid. We investigate the dynamics by measuring the self-intermediate scattering function 



Figure 6-a shows the evolution of F s (q,t) from the sol to the gel phase. While relaxation is fast and exponential in 
the sol phase, a slow secondary relaxation emerges at percolation. The final decay time varies little in the gel phase, 
but the height of the plateau at intermediate times changes strongly across percolation. A similar behavior is found 
for the coherent scattering function, as in experiments (13). Physically the plateau reflects the thermal vibrations 
of an elastic solid on intermediate timescale, while the long-time decay reflects the flow of the system: The system 
is viscoelastic. In Fig. 6-b we show that viscous flow is mostly controlled by rii n k, the rate for polymer extraction 
(4; 17). Flow in this system occurs when the percolating network slowly rearranges through polymer moves (14). 
Therefore, gelation corresponds to the continuous emergence, for increasing polymer concentration, of a plateau in 
dynamic functions, with an almost constant relaxation timescale, controlled by the polymer dynamics. Gelation is 
thus qualitatively different from a glass transition where the plateau height remains constant with a dramatic increase 
of relaxation timescales (18), though in both cases no structural signature of the transition is found. Coincidence 
of gelation and percolation, put forward in (4) or dispelled in (9), happens whenever long-lived bonds make cluster 
restructuration very slow, but does not occur in systems where the bond lifetime is short at percolation (6). Notice 
that the height of the plateau in F s (q,t) for small q as a function of R can be used as an order parameter for the 
gelation transition. 

Surprisingly, the mean-squared displacements shown in Fig. 6-c and 6-d appear as very poor indicators of the 
dynamics. The comparison with the behavior of the self-intermediate scattering function is in fact quite striking. 
While the final relaxation timescale, r, in F s (q,t) scales roughly as rii n k, the diffusivity, D s , extracted from the 
long-time behavior of A 2 (t) ~ QD s t is almost constant. This is reminiscent of the "decoupling" phenomenon, or 
"breakdown" of the Stokes-Einstein relation, reported in supercooled fluids (18). While "fractional" breakdown is 
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FIG. 6 Self-intermediate scattering function for q — 0.46 (a, b) and mean-squared displacement (c, d) for N = 10 3 , 4> = 0.2. 
(a, c) present the dynamics for rii n k = 10 3 and several values of R across percolation (R p « 1.85 for <f> = 0.2). Viscoelasticity 
continuously emerges at percolation, (b, d) are for R = 4 and different values of THnk, which directly controls the long-time 
decay of F s (q,t), while A 2 (t) remains essentially unchanged because its behavior is dominated by fast particles not belonging 
to the percolating cluster. 



reported in liquids, D s ~ , with £ in the range 0.7-0.9 instead of the normal value £ = 1 (19), we find here ( w 0, 
quite an extreme case of decoupling. Decoupling in gels was reported in different systems before (20). 

In supercooled fluids, decoupling phenomena are commonly attributed to the existence of dynamic heterogeneity, 
that is, the existence of non-trivial spatio-temporal distributions of mobilities. The analogy is confirmed in Fig. 7 
where we show distributions of droplet displacements, 

G.(r,t) = (± J2S(\r - ri (t) + Ti(0)\)^ . (3) 

Clearly, G s (r, t) exhibits a bimodal character suggesting coexistence of slow arrested droplets and fast diffusing 
droplets. Qualitatively similar distributions were reported in gels (7; 11; 21; 22) and glasses (19; 23; 24). Here, the 
snapshots in Figs. 3-e and 3-f suggest an obvious explanation for dynamic heterogeneity. At any given time, a system- 
spanning cluster of droplets which behaves as a solid on timescales smaller than Ti in ]< coexists with droplets which 
can diffuse through this arrested structure. We quantitatively confirm this interpretation in Fig. 7 where G g (r,t) is 
decomposed over two families of droplets: G s (r, t) — caGa(t, t) + (1 — ca)Gm(t, t), where A and M stand for droplets 
that are Arrested and Mobile at time t = 0, respectively; ca is the fraction of droplets belonging to the percolating 
cluster. While the central peak in G s (r,t) is dominated by GA(r,t), the large "non-Gaussian" tails are dominated by 
G M (r, t). 

A simple analytical model can be proposed to describe the dynamic heterogeneity of the equilibrium gel, which 
incorporates the physical idea of the coexistence between a slow, percolating cluster of connected droplets and fast, 
more freely diffusing droplets, with a dynamic exchange between the two families set by polymer moves. Similar 
physical ideas were qualitatively discussed earlier (11; 20; 25), but were, however, not exploited within a quantitative 
model. The proposed model, described in detailed in Ref. (3), can be solved analytically in the Fourier-Laplace domain 
and yields quantitative predictions for G s (r, t) in terms of a few free parameters (as for instance the concentration ca 
of arrested droplets, or the diffusivity of mobile particles) which have a clear physical interpretation and can be fixed 
by simple numerical observations. The lines in Fig. 7, which agree nicely with measured data, are predictions derived 
from this model. We thus find that the diffusion constant is in fact entirely dominated by those droplets which do not 
contribute to viscoelasticity, and is therefore a poor indicator of the gel dynamics. These results show that dynamic 
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FIG. 7 Distribution of droplet displacements for <j> = 0.2, R = 4, nink = 100 and different times for all droplets (O)i an< i its 
decomposition over droplets that are free (o) or arrested (□) at time t = 0. Lines are from a simple theoretical model based in 
this decomposition and introduced in Ref. (3). 



heterogeneity in gels can be stronger than in supercooled fluids, but its origin is also much simpler: The structure of 
the system is heterogeneous, Fig. 3, while no such static structure needs to be invoked in glasses to account for single 
particle dynamic heterogeneity (24). 



V. CONCLUSIONS 

In this work we studied a simplified model of a specific material : a microemulsion of oil droplets in water mixed 
with telechelic polymers which can form bridges between the droplets. Although motivated by this material, the 
model for reversible gelation described here sheds light on the microscopic aspects of gelation and the heterogeneous 
dynamics of gel-forming systems. Our modelling approach neglects the solvent, and uses droplets and polymers as 
elementary objects, with a coarse-grained description of polymers as simple links inducing an effective attraction 
between the droplets they connect. This system can be efficiently simulated using a hybrid molecular dynamics / 
Monte Carlo scheme. 

The model exhibits different phases as the control parameters, <j> and R, are varied. We find in particular a sol phase 
with liquid-like structure and dynamics, a coexistence region with both nucleation and spinodal decomposition regimes 
ending at a liquid-gas phase transition, a gelation line determined by geometric percolation, and an equilibrium gel 
phase unrelated to phase separation. The global structure of the equilibrium gel is homogeneous and no structural 
differences with the sol phase are detected. Despite this apparent homogeneity, the stress is supported by a fractal 
network which endows the system with viscoelastic properties. Moreover, cluster statistics and geometry at gelation 
are compatible with the random bond percolation universality class. 

We also show that, as opposed to structure, dynamics is strongly affected by gelation, where a slow secondary 
relaxation appears in time correlation functions related with the polymer rearrangement timescale, fUnk- I n this way, 
dynamical (but no structural) changes can be used to locate the gelation transition. The equilibrium gel dynamics 
is highly heterogeneous as a result of the presence of particle families with very different mobilities. In fact, at any 
given time, a system-spanning cluster of droplets which behaves as a solid on timescales smaller than rii n k coexists 
with droplets which can more freely diffuse through this arrested structure. These observation motivated a simple 
yet accurate analytic modelling of dynamic heterogeneity, which is generally applicable to gels (3). Although slow 
and heterogeneous dynamics are reminiscent of the physics of supercooled fluids, we discussed several qualitative 
differences between gels and glasses. 

In this paper we have discussed how gelation plays out at low volume fractions. At high volume fraction, crowding 
will result in glassy dynamics, and interesting questions arise when gel dynamics competes with glassiness. We are 
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currently investigating this issue. 
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